Effect of electron-electron interactions on the conductivity of clean graphene 
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Minimal conductivity of a single undoped graphene layer is known to be of the order of the 
conductance quantum, independent of the electron velocity. We show that this universality does 
not survive electron-electron interaction which results in the non-trivial frequency dependence. We 
begin with analyzing the perturbation theory in the interaction parameter g for the electron self- 
energy and observe the failure of the random-phase approximation. The optical conductivity is 
then derived from the quantum kinetic equation and the exact result is obtained in the limit when 
g <l<«?|lnw|. 

PACS numbers: 73.23.-b, 72.30.+q 



Introduction. Recent experiments on transport in 
graphene layers [UBS] have validated extensive theo- 
retical efforts directed at understanding of various prop- 
erties of two-dimensional Dirac fermions. A lot of these 
efforts are devoted to the zero temperature dc conduc- 
tivity which has a universal value of the order of the 
conductance quantum 0, [B}. Such a value is not un- 
expected from the dimension analysis since the intrin- 
sic (undoped) graphene lacks a characteristic momentum 
scale. Two different minimal conductivities are conven- 
tionally defined. The dc limit of an ac conductivity in 
a clean graphene (r _1 = 0, lu — > 0) was shown to be 
a = e 2 /Ah [6j. Another possible definition of a strict dc 



limit of impure graphene (lu = 0, r 
ent, but numerically close value a 



2e- 



0) gives a differ- 
1. Recent 



calculations have largely confirmed 12, IH 

the results of Ref. @, while others obtained different val- 
ues 1J] . Conductivity in a bilayer graphene has also been 
a subject of close theoretical attention [l7], 0, 0, 0|. 

It is noteworthy that the minimal conductivity a is 
very much analogous to the optical conductivity of a two- 
dimensional electron system with spin-orbit-split band 
structure, where it is due to the "chiral resonance" , and 
the corresponding value e 2 /16fi. [TBI, [l6| is exactly 4 times 
smaller than a (which is the degree of spin-nodal degener- 
acy in graphene) . This analogy is due to similarity in the 
chiral properties of the eigenstates in the two systems. 

In the present paper we demonstrate that the notable 
universality of the values of a and a does not hold in the 
presence of electron-electron interactions, which result in 
a strong frequency dependence of the conductivity. Here, 
we concentrate on the optical conductivity u(lu) in the 
strict disorder- free graphene (t _1 = 0) and show that the 
optical conductivity is actually suppressedby interactions 
in comparison with its "universal" value. 

Single intrinsic 2D graphene layer is described by the 
chiral Hamiltonian, 
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represents the usual set of Pauli matrices), the sum 
over Latin indices is taken over two nodal points and 
two (true) spin directions. The interaction potential is 
Vq = 27re 2 /K|q|, with k being the dielectric constant of a 
substrate, Vq is the "bare" electron velocity. Hereinafter 
we denote, J2 P = I d 2 p/(2n) 2 , and set % = 1 throughout 
the text, except in the final result ([20]) . 

Let us begin with the perturbation theory, in powers 
of the dimensionless interaction constant g — e 2 / kvq, for 
the electron self-energy (at T = 0). We find that the 
random phase approximation (RPA) fails for the system 
described by the Hamiltonian (TT]), as the non-RPA con- 
tributions are generally not small. 

Perturbation theory. The first and second order inter- 
action corrections to the electron self-energy are shown in 
Fig. 1. The solid line corresponds to the electron Green's 
function 
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where <7 p = <x • n is the projection of the pseudo-spin 
operator onto the direction of the electron momentum 
n = p/p. The first-order contribution (Fig. la) is inde- 
pendent of the energy variable, 



p' 



^ p _p,GV = i <7p^y p _p,(n- n'). 
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The integral here diverges logarithmically at p' 3> p. 
This divergence is cut-off at the inverse lattice spacing, 
/C, leading to the logarithmic correction, 
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p <7 P — ln(fC/p). 
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where "hats" denote operators in pseudo-spin space (a 



Such logarithmic renormalization of the electron velocity 
was first discussed in Ref. [2l| within the scope of the 
renormalization group approach. 

The second order correction consists of three contribu- 
tions, Fig. lb, c, and d, of which the last one is identically 
zero. This fact can be understood upon noticing that the 
diagram Id is responsible for the renormalization of the 
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,<z> 



c) 



d) 



The leading logarithmic contribution into this expression 
comes from large p',p" 3> p. To this leading order, it is 
sufficient to expand the integrand to the linear power in 
p. As the relevant terms appear from the interaction 
potentials as well as from the expansion of both the nu- 
merator and denominator, the calculations are too cum- 
bersome to be presented here. Finally, 



FIG. 1: First (a) and second (b,c,d) order corrections to the 
electron self-energy. Solid line represents the electron Green's 
function ((2]), dashed line stands for the bare interaction po- 
tential Kj- Closed loop implies summation over nodal and 
spin degeneracy. 



Fermi level 22J and must be absent in clean graphene due 
to its electron-hole symmetry. The other two graphs give 
non-zero contributions. The RPA-type diagram Fig. lb 
gives 
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where the coefficient 4 comes from the summation over 
the nodal points and (real) spin states in the internal 
loop which represents the polarization operator 



n(w,g)=-iTr^ J ^Gv +W:P / +q G £V , 
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here trace is taken in the pseudo-spin space. To evaluate 
the integrals it is convenient to use the following integral 
representation for the polarization operator, 
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The integral over frequency in Eq. ([5]) can now be easily 
calculated. Restricting for simplicity to the e = case 
here 23JI, we obtain, 



dx 



p q 



P kvq^Jq y/l + x 2 qy/l + x 2 + |p - q| 



(8) 



The main (logarithmic) contribution to this integral 
comes from large values of transferred momentum q. Ex- 
panding to the linear order in p/q, and evaluating hrst 
the integral over the angle between p and q, then the 
integral over dx and finally over the absolute value of q, 
we arrive at 

c 4 



p ^ " p 6kvq 



ln(/C/p). 
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The remaining contribution Fig. lc is given by (again, 
e = 0), 
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The expressions (|4I9|11|) determine the remormalization 
of the electron velocity up to the second order in the 
electron-electron interaction, 
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Two important conclusions can be drawn from this result. 

(i) The contribution (fTTjl from the non-RPA diagram 
with intersecting interaction lines is not parametrically 
small compared with the RPA term ([5]). Resulting in 
an overestimation of about ~ 20% in the second order, 
the neglect of non-RPA corrections to the electron self- 
energy becomes an uncontrollable approximation in the 
higher orders in g for an undoped graphene. Our findings, 
thus, do not support the conjecture of Refs. (where 
it was used for the calculation of the electron lifetime) 
that RPA is exact approximation in the limit j C 1. 
Experimentally, the value of the velocity is vq = 1.1 x 10 6 
m/s [2j. For a typical dielectric substrate k rs 6, which 
yields g w 0.3, indicating that g <C 1 is a reasonable 
approximation. 

(ii) It is essential that the higher orders do not produce 
higher powers of hh/C/p). This is most simply verified 
for the RPA-terms [2l|. On the other hand, the diagrams 
with intersecting interaction lines are not more singular 
than the RPA ones as is clear from the power counting 
(but also not parametrically small than RPA terms) . As 
the consequence, in the limit g <C 1 it is sufficient to 
restrict to the lowest order of the perturbation theory. 
Still, the product g\n(JC/p) can be arbitrarily large with- 
out violating the perturbation expansion. 

Kinetic equation. We now apply the obtained under- 
standing to the analysis of the ac conductivity o~(u>). 
Free-electron conductivity follows directly from the po- 
larization operator {7J and the particle conservation con- 
dition. Taking into account nodal and spin degeneracy, 

■ 2 2 

%e uj 6 
<r ( w )=41im L_n(w,g) = T . (13) 
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In order to calculate the homogeneous ac conductivity 
in the presence of electron-electron interactions we apply 
a quantum kinetic equation. This will allow us to take 
into account the leading terms ~ g\n(JC/p) in the exact 
way (but still assuming g <C 1). The density matrix f p 
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should be defined as a matrix in a pseudo-spin space, 
fp 13 = (( c p tc p))j where a and (3 are the pseudo-spin in- 
dices. Using the equations of motion for the operators c p 
determined by the Hamiltonian ([1]), and assuming a ho- 
mogeneous in space (but time-dependent) external elec- 
tric field E, it is straightforward to write a closed equa- 
tion for the density matrix / p to the linear order in the 
electron-electron interactions, 



dfp 

dt 



+ ivp[a p ,f p ] + eE 



dp 



WLfp'Jp]- (14) 
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We emphasize that kinetic equation for graphene ([T] 
does contain the terms linear in the interaction. This 
is in contrast to the conventional system with parabolic 
dispersion for which non-trivial contributions (collision 
integral) appear in the second order in Vq [25]. This 
difference is due to the chiral structure of the electron 
eigenstates in graphene. Kinetic equation needs to be 
solved to the linear order in the applied electric field, 
f p = ,/p°^ + f p , around the equilibrium distribution 
fp ^ = 1 — <Tp/2 (at zero temperature). The contribution 
resulting from the term [f p /,fp] represents simply the 
electron self-energy to the first order in interaction. 

The remaining term [/!/ , /p ] describes vertex correc- 
tions. For the periodic electric field the linearized kinetic 
equation takes the form 

- iw/W + iv p p [dp, /£)] = ^(E x n) ■ {a x n) 



2p 
p' 



Vp-p'[£\a p ], (15) 



where the renormalized velocity v p is given by the expres- 
sion (fT2|) taken to the first order in g. With the following 
substitution, 

= i(E X n) • (<x X n)A(p) + E • (a x n)B(p), (16) 



the matrix integral equation (|T5|) reduces to a pair of 
coupled scalar integral equations, 

— ^ v p-p ,b (p') cos v 
p' 



;A(p) + 2v pP B(p) = 
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2v pP A(p)+uB(p) = ^Vp-p, A{p') cos 8 PP >, (17) 
p' 

where 9 pp i is the angle between the vectors p and p'. The 
system of equations (fl7|) can be solved by means of the 
consecutive approximations in Vp- P >. Fig. 2 illustrates 
diagrammatically the meaning of these consecutive or- 
ders. In the limit g <C 1 it is sufficient to stop at the first 
order term (see the discussion at the end of this section), 
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FIG. 2: Graphic representation of the density matrix given 
by the integral equation (|15[) . The dashed line represents 
the interaction potential F q , the solid line stands for the 
electron propagator corrected by the self-energy E™. The 
cross denotes the applied external field. The first order term 
brings a logarithmic contribution into the optical conductiv- 
ity, oc f dp/vpp. The higher-order terms contain higher power 
of the small coupling constant g but do not result in higher 
powers of logarithms. 



Here, as usual, the frequency should be understood as 
posessing the infinitesimally small positive imaginary 
part, uj — > lu + ir\. Knowledge of the function A(j>) allows 
to determine the electric current, j = 4eTr ^ p uojpd - = 
cr(w)E, which yields the conductivity 
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The real part of the conductivity is due to the zeros of 
the denominators of Eq. (fl8|) . In the second term in 
Eq. (|18p the leading logarithmic contribution comes from 
the imaginary part that arises due to the second singu- 
larity — 2v P 'p'). It is important that the second 
order term in the expansion Eq. (|18p . which is not ex- 
plicitly written, does not yield any g 2 hi 2 p contribution, 
being only ~ g 2 hip, and, thus, by virtue of g -C 1 result- 
ing only in a small correction to the first term, similar 
to the situation already encountered in the self-energy, 
cf. Eq. (|12p . Upon extracting this imaginary part, one 
should first evaluate the integral over the angle 9 pp i and 
then over the momentum p. Straightforward calculation 
leads to the low-frequency optical conductivity 
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where = Vo[l+gbx [jj^j /4] is the renormalized elec- 
tron velocity, corresponding to frequency to. In deriving 
Eq. ([20]) it is utilized that J^ /vo dp/v p p = 4 In ( %A /gv . 

Eq. ([20|) is the main result of this paper. The conduc- 
tivity depends only on the combination of the frequency 

and the coupling constant, gin (^jfjfj ■ We illustrate the 
dependence of the optical conductivity on the frequency 
for different values of the coupling constant g in Fig. 3. 
When frequency increases the optical conductivity ap- 
proaches its value e 2 /Ah, but at lower frequencies its 
magnitude is significantly reduced. 

Restricting the perturbation expansion ([T8|) to the low- 
est order vertex correction requires some justification. 
In particular, one has to make sure that this expansion 
does not lead to higher-order logarithms in the discarded 
terms, i.e. that no x-ray (Mahan) singularity [22j takes 
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FIG. 3: (Color online) The dependence (|20|) of the real part 
of the optical conductivity (measured in units of the non- 
interacting conductivity e 2 /4?i on the dimensionless frequency 
= uj/vqK, for different values of the coupling constant: 
g = 0.5 (green dotted line), g — 0.7 (red dashed line), g = 1.0 
(black solid line) . The suppression of the conductivity is pro- 
nounced in the low-frequency domain. 



place. Indeed, the second order contribution into A(p) 
can be easily derived and analyzed. The corresponding 
expression is cumbersome but the conclusion is straight- 
forward. The higher order terms still feature single- 
logarithmic divergencies but arrive with higher powers 
of g. The situation here is similar to the one already 
encountered above in the calculation of the electron self- 
energy, see Eq. (fT2|) . 

Summary and conclusions. Optical conductivity in 
undoped graphene is the result of the "crural" resonance, 
i.e. resonant creation of an electron-hole pair. Chiral 
structure of the Hamiltonian leads to non-zero matrix 
elements of the velocity operator for the inter-subband 
transitions, resulting in a finite conductivity in a homo- 
geneous electric field (2o| . 

Effects of electron-electron interaction are two-fold. 
First, the logarithmic renormalization of the electron ve- 
locity at low energies leads to the decrease in the den- 
sity of states, which suppresses the probability of inter- 
subband transitions. This effect is accounted for by the 
electron self-energy and reveals in the appearance of the 
factor Vq/v^ in the expression ([20")) . Second, the excited 
electron and hole interact in the final state. This inter- 
action is attractive and results in the relative enhance- 
ment of the optical conductivity. Such final-state inter- 
actions are analogous to the excitonic effect in conven- 
tional semiconductors and are accounted for by the ver- 
tex corrections, Fig. 2, yielding the factor 1+ln into 

Eq. (|20p . Utilization of the quantum kinetic equations al- 
lows to consider both these effects on equal footing and 
account for the leading logarithmic terms in the limit 
when g < 1 <C g In \^jf) ■ Note that in order to be able 
to study homogeneous conductivity one should be able to 
neglect charge accumulation at the boundaries and asso- 
ciated with it electric field. It is easy to estimate that 
the system size should be larger enough, L 3> gvo/uj. 
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